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ABSTRACT 


A  program  has  been  written  for  the  HP-97  (HP-67)  minicomputer  to  solve 
the  blast  wave  from  a  nuclear  detonation.  The  program  first  determines 
standard-altitude  parameters  for  the  altitude  of  interest.  Using  these 
parameters  as  input,  the  program  then  calculates  pressures,  temperatures, 
densities,  velocities,  and  Mach  Numbers  in  front  and  behind  the  shock. 

C<  07  = 


Three  subroutines  allows  the  user  to  input  (for  any  aititudc^^82,000  ft)  a 
specific  overpressure  across  the  shock,  a  specific  gust  velocity  behind  the 
blast  wave,  or  a  specific  dynamic  pressure  behind  the  blast  wave.  Corre¬ 
sponding  parameters  for  each  input  are  calculated. 
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NUCLEAR  BLAST  PROGRAM 


For  Mini-Calculators 

Introduction 

In  the  process  of  performing  calculations  of  the  values  of  the  pres¬ 
sures,  temperatures,  velocities,  densities,  etc.  behind  weak  shock  waves  at 
some  distance  fn  m  a  nuclear  detonation,  it  became  obvious  that  hand  calcu¬ 
lation  was  tedious  and  interpolation  from  graphical  results  subject  to 
e-ror.  Therefore  a  simple  program  was  developed  for  use  on  the  HP-97  (or 
HP-i7)  mini  calculator  to  perform  these-  operations.  !his  program  is 
described  in  the  following  narrative. 

Although  this  simple  program  is  not  intended  to  replace  the  more 
sophist icated  computer  codes  available  for  such  computation,  it  is  a 
valuable  tool  for  first  cut  analyses.  It  also  has  the  distinct  advantage  of 
being  highly  portable,  and  serves  as  handy  reference  for  use  when  the 
detailed  data  are  not  immediately  available. 

Discussion 

At  distances  from  any  nuclear  detonation  of  interest  to  the  analyst 
working  with  airborne  systems,  the  blast  wave  is  relatively  weak  and 
amenable  to  solution  via  standard  fluid  mechanical  equations.  The  first 
step  in  such  a  solution  is  tVie  application  of  standard  fluid  dynamical 
equations  used  to  solve  the  stationary  shock.  See  Figure  1.  The  governing 
equations  relating  conditions  in  front  of  and  behind  the  shock  can  be 
obtained  from  any  text  on  supersonic  flow,  I  used  Shapiro*. 


*  A.  H.  Shapiro,  The  Dynamics  and  Thermodynamics  of  Compressible  Fluid  Flow. 
Ronald  Pres®,  N.Y.  1953 
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Figure  1 . Flow  parameters  Cor  a 
stationary  shock  wave 


After  solving  the  stationary  shock  problem,  I  then  transformed  the 
problem  into  a  moving  shock,  or  blast  wave  probLem  by  allowing  the  shock  to 
move  with  a  speed  such  that  its  Mach  number  was  equal  to  Mj .  This  system 
uses  primed  variables  and  is  depicted  in  Figure  2. 
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I  Figure  2.  Flow  parameters  for  moving  shock  in 
terms  of  stationary  shock  solution. 


As  depicted  in  Figure  2,  the  stationary  shock  solution  provides  all  the 
data  needed  to  solve  the  blast  wave  problem.  The  parameters  of  moat 
interest  to  the  nuclear  survivability/vulneiability  analyst  are  Vj*  (the 
gust  velocity),  the  overpressure,  ( P2  *  — ■  Pj  ■  P2  ~  Pj )  and  the 
dynamic  pressure,  . 
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The  program  i ncorporat  wig  the  eolation  ot  ttu*  blast  wave  is  listed  in 
the  Appendix  H  and  detailed  operating  instructions  are  contained  in  Appendix 
A.  The  ''gut a”  of  the  program  is  subroutine  B,  which  solves  the  stationary 
shock  then  utilizes  that  solution  to  solve  the  blast  wave  (moving  shock)  for 
any  selected  altitude  and  overpressure.  Subroutines  C  and  D  accept  gust 
velocity  and  dynamic  pressure  inputs  respectively  and  iterate  back  through 
the  basic  part  o!  subroutine  B  until  the  calculated  quantity  is  within  the 
acceptable  tolerance  of  the  input. 

Subroutine  A  is  the  cornerstone  of  the  program.  Tt  contains  analytical 
expressions  describing  the  variation  of  the  standard  atmosphere  with 
altitude  as  presented  in  Dommasch.*  For  the  specified  altitude  these 
equations  are  solved,  resulting  in  ratios  of  the  values  of  the  parameter  of 
interest  at  altitude  to  the  sec  level  standard  value.  Multiplication  of 
these  ratios  by  the  pertinent,  sea  level  standard  results  in  the  desired 
value  at  altitude.  These  values  are  stored  and  also  used  in  subsequent 
calculations.  Note  that  because  of  program  step  limitations,  Che  program  was 
separated  into  two  parts.  Fact  1  ta  pertinent,  for  the  troposphere  (b< 

36,000  feet)  whi  le  part  2  i. ;  accurst.'-  in  the  stratosphere  (36,000<h  <_ 

82,000)  where  the  temperature  is  constant  with  altitude. 


CONCLUSION 

The  program  presented  herein  should  prove  to  be  a  valuable  tool  to  analysts 
dealing  with  weak  blast  waves  of  any  type,  but  should  prove  of  particular 
interest  to  nuclear  survivability/vulnerability  analysts  in  government  and 
industry. 


*D.  0.  Dommasch,  S.S.  Sherby,  and  T.F.  Connolly,  Airplane  Aerodynaai< 
Pitman  Publishing  Corp.  N.Y.  1961. 
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1.  Subroutine  A,  Calculation  of  Ambient  Atmospheric  Parameters. 

Upon  loading  the  program,  the  first  step  must  be  the  utilization  of 
Subroutine  A  to  calculate  the  ambient  static  pressure  and  temperature  at  the 
altitude  of  interest.  The  other  subroutines  use  these  parameters,  and  will 
not  function  unless  the  calculations  are  made.  Note:  Load  Part  1  if 
altitude  of  interest  is  36,000  feet  or  less,  load  Part  2  if  the  altitude 
of  interest  is  between  36,000  and  82,000  feet. 

PRESS  A  -  This  step  clears  all  registers  and  sets  up  the  program. 

ENTER  ALTITUDE  -  Altitude  is  in  feet. 

PRESS  R/S  -  Calculation  requires  a  few  seconds.  When  the  altitude  is 
displayed,  calculations  are  complete.  Secondary  registers  contain  results 
of  the  calculations. 
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To  recaLL  these  par  am-  ters,  press  t,  then  P**S.  Use  RCL  for  register  of 
Interest.  P»\  ->r  to  use  of  other  subroutines  either  restore  the  original 
configuration  ( 5  ,  P^S)  or  restart  program. 

II.  Subroutine  B  -  Calculation  of  Gust  Velocity  (given  an  Overpressure). 
This  subroutine  accepts  an  overpressure  (in  psi)  and  solves  the  blast 

wave.  The  gusi  velocity  is  displayed  at  the  completion  of  the  calculations. 
(Recall  that  Subroutine  A  must  be  exercised  at  least  once  prior  to  use  of 
this  subroutine.  The  results  of  the  subroutine  are  pertinent  only  for  the 
last  altitude  entered  in  Subroutine  A.) 

PRESS  R  Initializes  subroutine 

ENTER  OVERPRESSURE  in  psi.  Recall  overpressure  *  P2  *  P]^. 

PRESS  R/S  Alter  calculations  are  complete,  gust  velocity  is 

d  i.  apt  ayed . 

At  the  point  the  user  has  several  options.  If  the  gust  velocity  corre- 
spending  to  the  selected  overpressure  for  a  different  altitude  is  desired, 
press  A,  enter  the  now  altitude,  press  R/S,  then  exercise  subroutine  B 
again,  Tf  uac ious  oveivressiu  e/gust  velocity  pairs  are  desired  for  the 
selected  altitude,  it  is  not  necessary  to  recalculate  the  atmospheric  param¬ 
eters,  press  B,  then  enter  the  overpressure.  The  corresponding  guat  is 
displayed . 

III.  Subroutine  C.  Calculation  of  Overpreiaure  (given  a  guat  velocity). 
This  subroutine  accepts  a  gust  velocity  (in  feet/second)  and  calculates 

the  corresponding  overpressure  via  iteration  through  the  basic  solution 
incorporated  in  subroutine  B.  This  subroutine  may  be  exercised  immediately 
after  the  altitude  computations,  or  the  completion  ot  any  of  the  other 


subroutines. 


PRESS  C 


Initializes  subroutine 


ENTER  GUST  VELOCITY  feet/second 

PRESS  R/S  Several  iterations  may  be  necessary  for 

convergence.  After  the  calculations  are  complete, 
the  overpressure  is  displayed. 

If  other  gust  velocity/overpressure  pairs  are  desired  for  the  same  altitude, 
simply  repeat  the  above  steps.  If  a  new  altitude  is  desired,  subroutine  A 
must  be  reaccomplished  prior  to  initialization  of  subroutine  C. 

IV.  Subroutine  D.  Calculation  of  the  Overpressure  and  Guat  Velocity 
(given  a  dynamic  pressure). 

This  subroutine  accepts  a  dynamic  pressure  in  psi  and  calculates  the 
corresponding  overpressure  and  gust  velocity  via  interation.  This 
subroutine  may  be  exercised  immediately  after  the  completion  of  subroutine 
A,  or  upon  completion  of  any  of  the  other  subroutines. 

PRESS  I)  Initializes  subroutine, 

ENTER  DYNAMIC  PRESSURE  in  psi.  NOTE.  A  maximum  of  3  places  after  the 

decimal  will  be  tolerated  -  otherwise  program 

cycles  endlessly. 

PRESS  R/S  Several  iterations  may  be  required  for 

convergence  after  calculations  are  complete, 
overpressure  is  displayed, 

PRESS  R/S  Gust  velocity  is  displayed. 

V.  Miscellaneous  Notes 

1.  If  the  program  cycles  endltssly,  or  if  an  error  is  indicated  in  the 
display,  press  A  twice  to  clear. 
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2.  The  ok* rci  r.e  t>f  •3'tl-ro'.,r  >  nos  B , 
l^^ne re ted  >  I ti  each  c a r e  ,  only  one 
but  oU'-i: '■  rthccV  ..'tiuns  ait 

needed.  The  register  contents  are 
figures  3  and  2. 


iod  T)  results  in  much  data  being 
(or  two  for  D)  parameters  are  displayed,'1 
■itore*  and  may  be  recalled,  if  they  are 
depicted  below.  Variable  labels  are  from 
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